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' Abstract 



Whightman function, the vacuum expectation values of the field square and the energy- 
momentum tensor are investigated for a scalar field inside a wedge with and without a 
£S| '. coaxial cylindrical boundary. Dirichlet boundary conditions are assumed on the bounding 

surfaces. The vacuum energy-momentum tensor is evaluated in the general case of the 
curvature coupling parameter. Making use of a variant of the generalized Abel-Plana formula, 
expectation values are presented as the sum of two terms. The first one corresponds to the 
geometry without a cylindrical boundary and the second one is induced by the presence 
of this boundary. The asymptotic behaviour of the field square, vacuum energy density 
, and stresses near the boundaries are investigated. The additional vacuum forces acting on 

the wedge sides due the presence of the cylindrical boundary are evaluated and it is shown 
that these forces are attractive. As a limiting case, the geometry of two parallel plates 
perpendicularly intersected by a third one is analyzed. 

■ PACS numbers: ll.10.Kk, 03.70.+k 
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1 Introduction 

The influence of boundaries on the vacuum state of a quantum field leads to interesting physical 
consequences. The presence of reflecting boundaries alters the zero-point fluctuations spectrum 
and results in the shifts in the vacuum expectation values of quantities quadratic in the field, 
such as the energy density and stresses. In particular, vacuum forces arise acting on constraining 
boundaries. The particular features of these forces depend on the nature of the quantum field, 
the type of spacetime manifold, the boundary geometries and the specific boundary conditions 
imposed on the field. Since the original work by Casimir many theoretical and experimental 
works have been done on this problem (see, e.g., |H 01 EJ El and references therein). How- 
ever, there are still difficulties in both interpretation and renormalization of the Casimir effect. 
Moreover, the absence of a complete renormalization procedure, in practice, limits all exact 
calculations to the special case of highly-symmetric boundary configurations (parallel plates, 
sphere, cylinder) with a specific background metric. From this point of view the wedge with a 
cylindrical boundary is an interesting system, since the geometry is nontrivial and it includes 
two dynamical parameters, radius and angle, for phenomenological purposes. Due to the formal 
analogy that exists between a wedge and a straight cosmic string, the corresponding results can 
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be applied to cosmic strings. In this paper we will study the vacuum expectation values of the 
field square and the energy-momentum tensor for a massless scalar field with a general curvature 
coupling parameter inside a wedge of opening angle 4>q with and without a coaxial cylindrical 
boundary assuming the Dirichlet boundary conditions on the constraining surfaces. The both 
regions inside and outside the cylindrical shell are considered. In the discussion below we will 
call these regions as region I and region II (see figure^). Some most relevant investigations to the 
present paper are contained in El El E] for a conformally coupled scalar and electromagnetic 
fields in a four dimensional spacetime. The interaction between minimally and conformally cou- 
pled scalar fields and different boundary configurations constructed from parallel and orthogonal 
plane surfaces is studied in The total Casimir energy of a semi-circular infinite cylindrical 
shell with perfectly conducting walls is considered in JJ| by using the zeta function technique. 
For a scalar field with an arbitrary curvature coupling parameter the vacuum energy density 
in the geometry of a wedge with an arbitrary opening angle and with a cylindrical boundary 
is evaluated in ^2] • This geometry is interesting from the point of view of general analysis for 
surface divergences in the expectation values of local physical observables for boundaries with 
discontinuities. The nonsmoothness of the boundary generates additional contributions to the 
heat kernel coefficients (see, for instance, the discussion in |13l 1141 ITd] and references therein). 
As in ^2j, our method here employs the mode summation and is based on a variant of the 
generalized Abel-Plana formula ^S] combined with the point-splitting technique. This allows 
to extract from the vacuum expectation values the parts due to a wedge without a cylindrical 
shell and to present the cylindrical parts in terms of strongly convergent integrals. 

We have organized the paper as follows. The next section is devoted to the evaluation of 
the Wightman function for a massless scalar field with a general curvature coupling inside a 
wedge with a cylindrical boundary. By using the formula for the Wightman function, in section 
131 we evaluate the vacuum expectation values of the field square and the energy-momentum 
tensor inside a wedge without a cylindrical boundary. In section |3] we consider the vacuum 
densities for a wedge with the cylindrical shell. Formulae for the shell contributions are derived 
and the corresponding surface divergences are investigated. The Wightman function, vacuum 
expectation values of the field square and the energy-momentum tensor in the region II are 
considered in section El As a special case of the considered geometry, in section El we discuss the 
limiting case of two parallel plates perpendicularly intersected by the third plate. Finally, the 
results are summarized and discussed in section 

2 Wightman function inside a wedge 

To describe the geometry of a wedge with the opening angle cfto and with the cylindrical 
boundary of radius a (see figure we will use the cylindrical coordinates (x 1 ,^ 2 , . . . ,x D ) = 
(r, <f>, Zi, . . . , zn), N = D — 2, where D is the number of spatial dimensions. Consider a massless 
scalar field <p with the curvature coupling parameter £, satisfying the field equation 



We quote the generalization to other boundary conditions in section |7| In (|2.1[) V* is the 
covariant derivative operator, R is the scalar curvature for the background spacetime. The values 
of the curvature coupling parameter £ = and £ = (D — 1)/4D correspond to special cases of 
minimally and conformally coupled scalars respectively. In this section we evaluate the positive 
frequency Wightman function (0\<p(x)tp(x')\0) in the region I, with |0) being the amplitude for the 




(2.1) 



and obeying Dirichlet boundary condition on the bounding surfaces: 
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corresponding vacuum state. The response of a particle detector in an arbitrary state of motion 
is determined by this function (see, for instance, ^ZI)- I n addition, the vacuum expectation 
value of the energy-momentum tensor is expressed in terms of the Wightman function as 



(0\T lk (x)\0) = lim ViV k (0\tp(x)<p(x')\0) + 

X'— >x 

Here we have assumed that the background spacetime is flat and have omitted the term with the 
Ricci tensor. By expanding the field operator and using the standard commutation relations, 
the Wightman function is presented as the mode sum 

(0|99(x)^(x')|0) = <Pa(*)<P%W)> ( 2 - 4 ) 

a 

where {(p a (x)} is a complete orthonormal set of positive frequency solutions to the field equation, 
satisfying the corresponding boundary conditions, a is a set of corresponding quantum numbers. 




\ 



4 ) 9ikV l Vi 



<<V(s)|0>. (2.3) 



Figure 1: Geometry of a wedge with the opening angle </>o and cylindrical boundary of radius a. 

In the region < r < a (region I in figure the eigenfunctions satisfying the boundary 
conditions (|2.2|) on the wedge sides <fi = 0, 0o have the form 

<p a (x) = l3 a Jqn(ir)s'm(qn(j))exp (ikry - iujt) , a = (n,7,k), (2.5) 
uj = a/7 2 + k 2 , q = ir/<f)o, —oo<kj<oo, n = 1,2, (2.6) 

where k = (ki, . . . , fcjv), ry = (zi,...,zn), and Ji(z) is the Bessel function. The normalization 
coefficient f3 a is determined from the standard Klein-Gordon scalar product with the integration 
over the region inside the wedge and is equal to 



fJa (270^^ J^( 7 a)' V ' ; 

The eigenvalues for the quantum number 7 are quantized by the boundary condition (|2.2|) on 
the cylindrical surface r = a. From this condition it follows that the possible values of 7 are 
equal to 

l = \ n ,j/a, J = 1,2,- • • , (2.8) 

where A nj - are the positive zeros of the Bessel function, J qn (^n,j) = 0, arranged in ascending 
order, \ n> j < Xnj+i- Substituting the eigenfunctions (|2.5|) into mode sum formula (|2.4|) with the 
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set of quantum numbers a = (n,j, k), for the positive frequency Wightman function one finds 



(0\<p(x)<p(x')\0) 



d N \t e* k ^ r " r H ^ sin(qn(f)) sm(qruf)' 
n=l 

oo 

^^J (?n ( 7 r)J gn ( 7 r')e--( i -*') 



(2.9) 



This formula is not convenient for the further evaluation of the vacuum expectation values of 
the field square and the energy- momentum tensor. This is related to that we do not know the 
explicit expressions for the eigenvalues A n as functions on n and j, and the summands in the 
series over j are strongly oscillating functions for large values j. In addition, the expression 
on the right of ()2.9|) is divergent in the coincidence limit and some renormalization procedure 
is needed to extract finite result for the vacuum expectation values of the field square and the 
energy- momentum tensor. 

To obtain an alternative form for the Wightman function we will apply to the sum over j a 
variant of the generalized Abel-Plana summation formula 16 



2/(A Hi j 



1 ^n,jJ'qn(^n,j)\/^nj + C 
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Vz 2 + c 2 4 



f(z)Y qn (z) 
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7T 



Vz 2 + c 2 J qn {z) 
Kq n (z) e- qnm f(ze^) + e^/Cze-Tf) 







Iqn(z) 



v 7 ^ 



+- 

7T 



" dz K ^ z~ qnm f(^) ~ e^fjze-^) p 1Q) 



Iqn\z) 



Vz^ 



where Y[{z) is the Neumann function, and Ii(z), K[(z) are the Bessel modified functions. This 
formula is valid for functions f(z) analytic in the right half-plane of the complex variable z = 
x + iy and satisfying the conditions 



\f(z)\ <e(x)e c ^, Cl <2, 



(2.11) 



f(z) = o(z 2 ^- 1 ), z^O, 



(2.12) 



where e(x) — > for x — > oo. By taking in formula (|2. 1U|) qn = 1/2, as a particular case we 
receive the Abel-Plana formula (see, for instance, [H). The generalized Abel-Plana formula was 
applied previously to a number of Casimir problems for spherically ^Hj and cylindrically |191 12(J| 
symmetric boundaries, and in the braneworld scenarios |21j . 

To evaluate the sum over j in ()2.9|) as a function f(z) we choose 



f(z) = zJ qn (zr/a)J qn (zr' /a) exp -i^/k 2 + z 2 /a 2 (t - t') 



(2.13) 



Using the asymptotic formulae of the Bessel functions for large values of arguments when n is 
fixed (see, e.g., [22]), we can see that for the function f{z) from (|2.13|) the conditions 1)2. 11|) . 
(I2.12|) are satisfied if r + r' + \t — t'\ < 2a. In particular, this is the case in the coincidence 
limit t = t' for the region under consideration, r, r' < a. Formula (|2.10|) allows to present the 
Wightman function in the form 



(0|^(x)^(x')|0) = (O w \tp(x)(p(x')\Q w ) + ((p(x)ip(x')) a , 



(2.14) 
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where 



i r d N k ik(r|l _ rli) r , ze -i(t-t')v^ 



(O w \<p(x)cp(x')\O w ) = - / -^ e ^ r l|- r ll^ / dz 



<f> J {2tt) n J V^ 2 + k 2 

oo 

^ sin(gn</>) sin(qn(j)')J qn (zr)J qn (zr'), (2.15) 



oo 

X 
n=l 



and 



jN], . poo zcosh (t — t')\J z 2 — k 2 

a «k(rii— rj.) / 



^>o J (2tt) n J k ^z 2 - k 2 

00 K (za) 

x ^ sm(qn<f>) sm(qn(j)')I qn (zr)I qn (zr')—^- — -. (2.16) 

n=1 I qn(Za) 

In the limit a —* 00 for fixed r, r', the term {ip(x)(p(x')) a vanishes and, hence, the term 
w ) is the Wightman function for the wedge without a cylindrical boundary with 
the corresponding vacuum state \0 W ). Consequently, the term (ip(x)(p(x')} a is induced by the 
presence of the cylindrical boundary. Hence, the application of the generalized Abel-Plana for- 
mula allowed us to extract from the Wightman function the part due to the wedge without a 
cylindrical boundary. For the points away the cylindrical surface the additional part induced by 
this surface, formula ([2.16ft . is finite in the coincidence limit and the renormalization is needed 
only for the part coming from the term (|2.15|) . 

3 Vacuum expectation values inside a wedge without a cylin- 
drical boundary 

In this section we consider the geometry of a wedge without a cylindrical boundary. Expression 
()2.15|) for the corresponding Wightman function can be simplified by using the formula [23] 



cos nb J qn {u) J qn {v) = ^- a n J I J u 2 + v 2 - 2wcos 2?m + - J , (3.1) 

n=0 ^ n=n- \ V * / 

where n± = ± [(qir =p 6)/(27r)] (square brackets mean the integer value of the enclosed ex- 
pression), a n = 1 for n / n±, a n± = 1/2 for \(qir =p b)/{2i:)\ = 0,1,2,..., a n± = 1 for 
\(qir =F 6) / (27r) I 7^ 0, 1,2, . . ., and the prime means that the summand with n = should be 
taken with the weight 1/2. This formula is valid for integer values of q = vr/</>o ^ 1. The results 
for other values q are obtained from the formulae for the renormalized vacuum expectation val- 
ues derived below in this section by analytic continuation over c/)q. Making use formula (|3.1|) . 
for the sum over n in formula (|2.15j) one finds 

00 ^2 n j+ 

^2sm(qn(/))sm(qn(/)')J qn (zr)J qn (zr') = — ^(-l) j+1 ^ a n Jo(zuj), (3.2) 

n=l j=l n=rij- 

where 



n j± 



r 2 _|_ r /2 _ 2rr ' cos[2n0 o + 4> + (-1) J >'], (3.3) 
7TT^+(-l)V)" 



20o 



(3.4) 
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Now, evaluating the integral over z in (|2.15|) with the help of formula 

~i(t-t')y/z 2 +k' 2 



f 
Jo 



dz- 



ze 



kJul-it-V) 3 



-Jo{zuj) 



v^ 2 + k 2 ^ 2 -(t-t') 2 ' 
the expression of the Wightman function is presented in the form 



(3.5) 



(0 w \cp(x)ip(x')\0 w } = ^ 



i 2(2vr) 



AT+l 



V u ? - (* - *o s 



(3.6) 



For the integral in this formula one has 



2 7V 7r (7V-l)/2 r (JV±1) 



u) - (t - t') 2 



n 2 + |r,| -rj,! 2 - (t - t') 2 



(JV+l)/2 ' 



(3.7) 



with r(z) being the Euler gamma function. Hence, for the Wightman function in the geometry 
of a wedge without a cylindrical boundary one obtains the formula 



(0 w \p(x)(p(x')\0 v 



'D-1A 2 %'+ 

rE E 



(-i) j+ V 



A 2+ 



|« 2 + |r|| - r' ]{ \ 2 - (t - t') 2 



(D-D/2- 



j=ln=rij- ^ 

Note that for the Wightman function of the Minkowski space without boundaries one has 



(3.8) 



(0m\<p(x)<pW)\0m) 



4tt 2 [|r-r'| 2 - (t - t') 2 } 



A 2 i(D-l)/2 : 



(3.9) 



where \0m) is the amplitude for the vacuum state in the Minkowski spacetime without bound- 
aries. This function coincides with the j = 1, n = term in formula ()3.8|) . Taking the coincidence 
limit x' — > x, for the difference of the vacuum expectation values of the field square, 



we find 



(^ 2 W)S = (0 w \^ 2 {x)\0 w ) - (0 M \<p 2 {x)\0 M ), 

r(V) ^ ^ (-iy+V 



(47T) 2 r 



2 J \ ^ \ ^ 



sin 



<S3)\D-1* 



(3.10) 



(3.11) 



j=l n=rij- 

where the prime means that the term j = 1, n = has to be omitted, and we use the notations 



n0 n + 



l + (-lp 



7TT^(l + (-l) J ") 

2(f) 



(3.12) 
(3.13) 



In Eq. (jSHJ), a n = 1 for n ^ n j± , a^ ± = 1/2 for |vr T ^(1 + (— l)- 7 ") |/20 o = 0, 1, 2, . . ., c%. ± = 1 
for 1 7r =F </>(! + (— 1) J ) |/2</>o 7^ 0, 1, 2, . . .. In the case -D = 3, by using the formula 



m— 1 



E 2 / n7T\ o 1 I 77171" \ 

sec H j = m esc ^rax H — — J , 



(3.14) 



n=0 
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for the renormalized vacuum expectation value of the field square one finds [7\ 



<^ 2 (*)> 



M 

rcn 



q 2 - 1 - 3g 2 esc 2 (q<j>) 



48vr 2 r 2 



(3.15) 



m 



0,1 



0) the main contribution in (|3.11|) comes 



Near the wedge boundaries <fi = 
from the terms j = 2, n = and n = — 1 for m = and m = 1 respectively, and the renormalized 
vacuum expectation value of the field square diverges with the leading behaviour 



r(^i) 



(4vr 



t (r\4> — 



(3.16) 



The surface divergences in renormalized vacuum expectation values of the local physical observ- 
ables are well known in quantum field theory with boundaries and result from the idealization 
of the boundaries as perfectly smooth surfaces which are perfect reflectors at all frequencies. 
These divergences are investigated in detail for various types of fields and general shape of 
smooth boundary 0121]. Near the smooth boundary the leading divergence in the field square 
varies as (D — l)th power of the distance from the boundary. It seems plausible that such 
effects as surface roughness, or the microstructure of the boundary on small scales (the atomic 
nature of matter for the case of the electromagnetic field [22]) can introduce a physical cutoff 
needed to produce finite values of surface quantities. An alternative mechanism for introducing 
a cutoff which removes singular behavior on boundaries is to allow the position of the boundary 
to undergo quantum fluctuations [201 ■ Such fluctuations smear out the contribution of the high 
frequency modes without the need to introduce an explicit high frequency cutoff. 

Now we turn to the vacuum expectation values of the energy-momentum tensor. By making 
use formula ()2.3j) . for the non-zero components one obtains (no summation over i) 



•i\(w) 
i I ren 



r( 



D+V 



2 D+2 



~2~; \ ^ 



tt 2 r 



n j + 

E 

j=ln=nj- 



(i) 



sin 



lU)\d+i 



\- t 2 /ren 



— £c)r(^ri) ^ a n sin <fffl cos 4>n ^ 



2 D^ r D 



sin 



where i = 0, 1, . . . , D, and we use the following notations 

/g = 1 + (4£ -1)[(D- 1)6^ sin 2 tf>® + D6 j2 

/j° ) -4D(e-e c ) S m 2 #, 



d 2 )\D+i 



0,3, 



D, 



f (i) 
Jjn 

f (2) 
Jjn 



D 



4 sin 



(3.17) 
(3.18) 



(3.19) 
(3.20) 
(3.21) 



In the case 4>q = ir/2 and for minimally and conformally coupled scalar fields, it can be checked 
that from formulae (|3.17j) and (|3.18j) . after the transformation from cylindrical coordinates to 
the cartesian ones, as a special case we obtain the result derived in jlOj . For a scalar field with 
a general curvature coupling parameter the vacuum energy-momentum tensor is non-diagonal. 
For the conformally coupled scalar field this tensor is diagonal and does not depend on the 
angular coordinate <fi. From the continuity equation VfcTj fc = for the energy-momentum tensor 
one has the following relations for the components: 



d r (rT-j 1 ) + rd^Tl 
d r (rTl) + rc^T 2 



0. 



(3.22) 
(3.23) 
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In the case of a conformally coupled scalar field we have an additional zero trace condition. It 
can be checked that vacuum expectation values (|3.17f) . ()3.18j) satisfy equations (|3,22|) . (|3,23|) , 

For a non-conformally coupled field vacuum expectation values (|3.17|) diverge on the bound- 
aries (j) = <j) m . For the points near these boundaries the leading terms in the corresponding 
asymptotic expansions are given by relations (no summation over i, i = 0, 3, . . . , D) 



lrpi\{w) 

\ I /ren 



/rp2\(w) 

\J-2 /ren 



"\ J 2 /ren 



out - Cc)r m 



2 D ir^ (r\(j)-(p r , 



\D+1 



(3.24) 



For the points away from the edge r = 0, this leading behavior is the same as that for the 
geometry of a single plate. 

In the most important case D = 3, by using formula (|3.14f) . for the components of the 
renormalized energy-momentum tensor we find 



ren 



\ 3 / ren 



1 



32vr 2 r 4 1 45 



+ 3^ 



Q 2 ) (£ ~ &) 



-12 



ren 



1 



32vr 2 r 4 



(e-e c )g 2 

sin 2 (qcj)) 
45 



sin (qcfi) 



3 y 3 



(3.25) 



+ 12 



(e-^c)g 2 

sin 2 (qcj)) 



sin (qcj)) 3 



\ J 2 / ren 



\ J 2 /ren 



3(g-&) g 3 cos (gg) 
87r 2 r 3 sin 

^(e-ec 



8vr 2 r 



2„4 



60 



1 



+ 



3g 2 



sin ((70) 



(3.26) 
(3.27) 
(3.28) 



For a conformally coupled scalar field this tensor coincides with the result previously obtained 
in literature [HI Ej . The corresponding vacuum forces acting on the wedge sides are determined 
by the effective pressure — (T|) r ^. These forces are attractive for the wedge with q > 1 and are 
repulsive for q < 1. Formula (|3.25|) for the component (Tq)^ is derived in |I2j by a different way. 
For the electromagnetic field in D = 3 the vacuum expectation values of the energy-momentum 
tensor are considered in |H] for a perfectly conducting wedge and in [0] for a dielectric wedge. 
In both cases the vacuum energy-momentum tensor is uniform. This is a consequence of the 
conformal invariance of the electromagnetic field. 

4 Field square and the vacuum energy- momentum tensor in the 
region I 

We now turn to the geometry of a wedge with additional cylindrical boundary of radius a. Taking 
the coincidence limit x' — > x in formula (|2.14j) for the Wightman function and integrating over 
k with the help of the formula [20] 



d N k 



k s g(z)dz 
IP 



7T N / 2 ( N + S I 

T(N/2) B [~^~ , 2 



dz z 



N+s-l 



ff(z)> 



(4.1) 



where B(x,y) is the Euler beta function, the vacuum expectation value of the field square is 
presented as a sum of two terms 



((V|o> = (o w \v 2 K) + (^ 



(4.2) 
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where the part induced by the cylindrical boundary is given by formula 



n n 1--P 
2 3 - D 7T — 



'D-l 



sin 2 (qncj)) 



dz z 



D-2 



K qn {z) ■ 
-^gn (^) 



(zr/a). 



Note that this part vanishes at the wedge sides = m , < r < a. Near the edge r 
main contribution into {(p 2 ) a comes from the summand n = 1 and one has 

4r- 2 (g + l)sin 2 (#)^(g) f_r_\2« 

,2a/ 



(4vr)^r(^l) 



for r -C a, with the notation 



£>+2q-2 -^g( z ) 



(4.3) 
the 

(4.4) 
(4.5) 



The part {f 2 ) a diverges on the cylindrical surface r = a. Near this surface the main contri- 
bution into Q4.3JI comes from large values n. Introducing a new integration variable z — * nqz, 
replacing the Bessel modified functions by their uniform asymptotic expansions for large values 
of the order (see, for instance, HU), an d expanding over a — r, for \tp — 4> m \ 3> 1 — r/a to the 
leading order one finds 



(q/2a 



.D-l 



D+l 

7T— r 



2 / 



V D-1 



r/2 



n=l 



D-2 -2ng(l-r/a)Vl+z 2 



(47r) 2 (a 

This leading behavior is the same as that for a cylindrical surface of radius a. 

Similarly, the vacuum expectation value of the energy-momentum tensor for the situation 
when the cylindrical boundary is present is written in the form 

<0|T tt |0) = (O w \TikK) + (T ik )a, (4.7) 

where {Ti^) a is induced by the cylindrical boundary. This term is obtained from the correspond- 
ing part in the Wightman function, (ip(x)(p(x')) a , acting by the appropriate differential operator 
and taking the coincidence limit [see formula (|2.3|) ]. For the points away from the cylindrical sur- 
face this limit gives a finite result. For the corresponding components of the energy-momentum 
tensor one obtains (no summation over i) 



\D-1 



(4.6) 



(47T) 



D-l 
2 



dzz 



D 



K qn (z) 



{«S^[V(y)] " a { J n [I q n{y)] cos(2gn0)} , y = 



(4.8) 



2(4tt) 



D-l 
2 



qn sin(2gn0) 



71=1 



dz z 



D 



-l K qn(z) 
Iqn(z) 



Xlqn(y) 



^V(y) + (1-40/^(2/) 



with the notations 



,(±) 



[g{y)\ 



+ — 2 

f 



}' 



{^-l){g'\y)+g\y) 

g'\y) + -g(y)g'(y) - g\y) (l ± [l - 4^(1 t l)] 
y I 

(4£ - 1) [g' 2 (y) + g\y)] - %{y)g\y) + ^g\y) (4£ ± 1) 



y 



(4.9) 

(4.10) 
(4.11) 
(4.12) 
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for a given function g(y), i = 0, 3, ...,D. Formula (|4.8|) for the energy density (i = 0) was 
previously derived in ^j. In accordance with the problem symmetry, the expressions for the 
diagonal components are invariant under the replacement <j) — > <pQ — <j), and the off-diagonal 
component (T^a changes the sign under this replacement. Note that the latter vanishes on the 
wedge sides 4> = 4>m, < r < a and for <fi = <fio/2. On the wedge sides one has cos(2gn</>) = 1 and 
for the diagonal components of the energy-momentum tensor we obtain (no summation over i) 

" r fiff^l g [ dz ^W'^ (413) 

where A, = 4£ — 1, i = 0,1,3,... ,D, A<i = 1. In particular, the additional vacuum effective 
pressure in the direction perpendicular to the wedge sides, p a = — (T 2 2 ) a \(f>=4> m , does not depend 
on the curvature coupling parameter and is negative for all values < r < a. This means that 
the vacuum forces acting on the wedge sides = 4> m due to the presence of the cylindrical 
boundary are attractive. The corresponding vacuum stresses in the directions parallel to the 
wedge sides are isotropic and the energy density is negative for both minimally and conformally 
coupled scalars. It can be checked that expectation values ()4.8|) . (|4.9|) satisfy equations (|3,22|) . 
(|^.2«-{j) and, hence, the continuity equation for the energy-momentum tensor. Aiming to compare 
with the result for the energy-momentum tensor of a cylindrical shell with the radius a let us 
write down the corresponding formula, which is obtained from the general result of (2131 and has 
the form 

^ = 2{ nll%% S rdzzB^atfMzr/a)], (4.14) 



n=— oo 



with the same notations as in (|4.1()j) - (|4.12|) . 

For < r < a the cylindrical parts (|4.8j) and (|4.9j) are finite for all values < 4> < 0o, 
including the wedge sides. The divergences on these sides are included in the first term on the 
right-hand side of (|4.7|) corresponding to the case without cylindrical boundary. Near the edge 
r = the main contribution into the boundary part (|4.8jl comes from the summand with n = 1 
and one has 

q Lb} +) - BH cos(2#)l F D (q) 2(? _ 2 
(V)a « 1 -5+i — ,„ 7T ~ I Ur , (4-15) 



(4^)^a I? + 1 r(^i)r 2 (g) V2 ^ y 
/T u „ [g£ + (1 ~ 4Qg] sin(2#)F D (q) / r x2g-l 

(2)a ~ ^Vrf^ ' ( 6) 

with the notations 

B\ ±] = (4^ - 1)(1 ± 1), i = 0,3,...,Z>, (4.17) 
B[ ±] = 4£(l/g - 1 ± 1) T 1 + 1, = 4£(2 - 1/q) ± 1 - 1, 

and the function .Fd (</) is defined by formula (|4.5j) . 

The boundary part (Zf ) diverges on the cylindrical surface r = a. Introducing a new 
integration variable z — > nqz and by taking into account that near the surface r = a the main 
contribution comes from the large values of n we can replace the Bessel modified functions by 
their uniform asymptotic expansions for large values of the order. Expanding over a — r, to the 
leading order for the diagonal components one finds 

/•oo 

,(_) (z)cos(2gn<£)l . (4.18) 
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where 



= &M 1+ ihi[* ±1+ iD-wt-i) }}' (419) 

Af\y) = 1--^[j/ 2 ±1 + ^(1t1)], (4.20) 
1 ~\~ y 

4%) = 4^-l + T ^[y 2 (4C-l)+4C±l], (4.21) 
On the wedge sides one has cos{2qn<fr) = 1 and this yields 

where the coefficients Ai are defined in the paragraph after formula (|4.13|) . Summing over n, to 
the leading order over (a — r) _1 one finds 

A T (^-) 

(T!)aU^ m ~ - *) 2 ) — , r-a. (4.23) 
2(47r) 2 ( a - r) D+1 

It can be seen that for the non-diagonal component to the leading order one has (T 2 1 ) a ~ 

{ a - r y D . 

For the angles < <fr < (fro, by using the formula 

f;„^_^(z^^(_^_ 1 ), (4 , 4) 

n=l 

introducing a new integration variable y = 2-7r(l — r/a)Vl + z 2 /<fr$, and expanding over (1—r/a), 
one finds that for |</>— </> m | 3> 1 — r/a the leading contribution of the term with ^l^(^) dominates 
and to the leading order we find 



D(t - £ c )r 



■D+V 



2 njn + i)/{ a _ r) n + v i = 0,2, . . . ,D. (4.25) 

This leading divergence coincides with the corresponding one for a cylindrical surface of the 
radius a (see, for instance, \2U\). For the other components to the leading order one has (T^) a ~ 
(T 2 1 ) a ~ (a — r)~ D . In figures [21 and |3] we have plotted the dependences of the boundary- 
induced vacuum expectation values for the components of the energy-momentum tensor in the 
case of D = 3 conformally coupled scalar field as functions on coordinates x = (r /a) cos (fr, 
y = (r/a) sin <fr for a wedge with (fro = tt/2. In particular, the corresponding energy-density 
is negative everywhere and the vacuum forces determined by the component (T 2 2 ) a |^ = m are 
attractive. 



5 Whightman function, the field square and the vacuum energy- 
momentum tensor in the region II 

Now we turn to the region outside the cylindrical shell (region II in figure [J): r > a, < <fr < <f>Q. 
The corresponding eigenfunctions satisfying boundary conditions (|2.2|) are obtained from (|2.5j) 
by the replacement 

Jqniir) g qn (jr,-fa) = J qn (ir)Y qn (-fa) - J 9ri (70^(77-), (5.1) 
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Figure 2: Dependence of the energy density and radial stress for a D = 3 conformally coupled 
scalar field, induced by a cylindrical boundary of radius a, on coordinates r and <f> for the 
wedge with the opening angle (fio = ir/2. The variables on the axes are x = (r/ a) cos <p and 
y = (r/a) sin </>. 




Figure 3: The same as in figure for the off-diagonal 2-component and the azimuthal stress. 
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where Y qn (z) is the Neumann function. Now the spectrum for the quantum number 7 is con- 
tinuous. To determine the corresponding normalization coefficient f3 a , we note that as the 
normalization integral diverges in the limit 7 = 7', the main contribution into the integral over 
radial coordinate comes from the large values of r when the Bessel functions can be replaced by 
their asymptotics for large arguments. The resulting integral is taken elementary and for the 
normalization coefficient in the region r > a one finds 



Pi 



(2k) 



2-D 



7 



[Jf„(7a) +Yin(l a )] 



(5.2) 



Substituting the corresponding eigenfunctions into the mode sum formula Q2.4J1 . the positive 
frequency Whightman function is presented in the form 



<0|^(xMx')|0> = -j- 



0o J (2tt) 



d N k 
Y' 



ik(ri,-r' ) 



sm(qn(f)) sin(gn<//) 



71=1 



19qn(jr, ja)g qn (jr',ja) ex P 



V7 2 + k 2 



(5.3) 



To find the part in the Wightman function induced by the presence of the cylindrical shell we 
subtract from (|5,3|) the corresponding function for the wedge without a cylindrical shell, given 
by (|2.15|) . In order to evaluate the corresponding difference we use the relation 



g g n(7r,7a)g g n(7r / ,7«) j , w ( f) _ _ 1 UH^) , )H (s) ( f) 



(5.4) 



where Hqn(z), s = 1,2 are the Hankel functions. This allows to present the Wightman function 
in the form 1)2.14(1 with the cylindrical shell induced part 



(<p{x)<p{x')) a 



d N k 

20 O j (27T)"' 



ik(ri|-r| | ) 



sin((7n0) sin(</n<//) 



n=l 



E 



''•/7^^(7r)^(7r')- 



i(t'-i)Vy+F 



rr(s)/ \ ? n 



(5.5) 



On the complex plane 7 we can rotate the integration contour on the right of this formula by 
the angle k/2 for s = 1 and by the angle —it/2 for s = 2. The integrals over the segments (0, ik) 
and (0, —ik) cancel out and after introducing the Bessel modified functions we obtain 



(<p(x)ip(x')) a 



d N k 



ik(rn-r;) 



00 z cosh 
dz 



(t - t')Vz 2 - k 2 



7T0O J (2k) N 

00 

sin(gn0) sm(qn(j)')K qn (zr)K, 



qn 



zr 



n=l 



\J z 2 — k 2 

f) Iqn(za) 

K qn (za)' 



(5.6) 



Taking the coincidence limit for the arguments, from this formula we obtain the vacuum expec- 
tation value of the field square in the region r > a: 



(A 



)3-£> 



7T 2 



sin 2 (gn</>) 



dzz°- 2V(z) 



K gn {Z 



■K 2 n (zr/a). 



(5.7) 
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Comparing with formulae (|2.16|) and (|4.3|) , we see that the parts induced by the cylindrical shell 
in the regions I and II are obtained from each other by the interchange I qn (z) ^ K qn (z). 

As for the region I, vacuum expectation value (|5.7|) diverges on the cylindrical surface. The 
leading term in the corresponding asymptotic expansion near this surface is obtained from that 
for the region I, formula ()4.6|) . replacing (a — r) by (r — a). For large distances from the 
cylindrical surface, r > a, we introduce in (|5,7|) a new integration variable y = zr /a and expand 
the integrand over a/r. The main contribution comes from the n = 1 term. By taking into 
account the value for the standard integral involving the square of the MacDonald function, to 
the leading order one finds 

k +*) re +t) n °-^ 



7T— a D~lT (D-1 + 2q) T 2 (q) 



For the part in the vacuum energy-momentum tensor induced by the cylindrical surface in 
the region r > a, from l|2..3|) . (|5.6|) . (|5,7|) one has the following formulae 



D-l 



x {affifefo)] " a U^n(y)] cos(2gn0)} , y = ^, (5.9) 
= r 2 ffi"T. f>sin(2gnfl f°° dz z D ^ ^ 
xK gn (y) 



K qn (z) 



2 A Kqn {y) + (l-^)K' qn {y) 



(5.10) 



with the functions a^ n [g(y)} defined by (IQl - dlH^I . By the way similar to that used above 
for the vacuum expectation value of the field square, it can be seen that for large distances 
from the cylindrical surface, r 3> a, the main contribution comes from the term with n = 1 
and the components of the induced energy-momentum tensor behave as (T-) a ~ (a/r) D+1+2q , 
(T 2 1 ) a ~ (a/r) D+2q . Near the cylindrical surface the leading terms of the asymptotic expansions 
for the components of the energy-momentum tensor in the exterior region are obtained from 
the corresponding formulae for the interior region, formulae (|4.2(-i|) . 1)4.25(1 . replacing (a — r) by 
(r — a). As for the interior region, the vacuum forces acting on the wedge sides due to the 
presence of the cylindrical shell are attractive and the corresponding energy density is negative 
for both minimally and conformally coupled scalars. 



6 Limiting case 

In this section we consider a limiting case of previous results, obtained in the limit <f>Q — > 0, 
r, a — > oo, assuming that a — r and a<f>Q = b are fixed. This limit corresponds to the geometry of 
two parallel plates separated by a distance b, perpendicularly intersected by the third plate (see 
figure 0J). The vacuum expectation values of the energy-momentum tensor for this geometry 
of boundaries are investigated in for special cases of minimally and conformally coupled 
massless scalar fields. Formulae presented below in this section generalize this results for the 
case of an arbitrary coupling and give an alternative representation of the vacuum energy- 
momentum tensor. We introduce rectangular coordinates (x n ,x' 2 , . . . , x' D ) = (x, y, z\, ■ ■ ■ , z^) 
with the relations x = a — r, y = rcf) in the limit under consideration. Below the components 
of the tensors in these coordinates we will denote by primes. In this limit, from the quantities 
corresponding to the wedge without a cylindrical surface we obtain the vacuum densities in the 
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region between two parallel planes. These quantities are well-investigated in literature (see [27] 
for the case of general mixed boundary conditions) and below we will consider the additional 
part induced by the presence of the intersecting plate at x = 0. The corresponding vacuum 
expectation values are obtained from the expectation values (• • ■ } a investigated in previous 
section. For this we note that in the limit under consideration one has q = vt/c/jq — > oo, and the 
order of the Bessel modified functions in formulae (|4.3|) . 1)4.8(1 . ((4.9)1 tends to infinity. Introducing 
a new integration variable z — > qnz, we can replace these functions by their uniform asymptotic 
expansions for large values of the order. The vacuum expectation value of the field square is 
presented in the form 

(0| ( , 2 |0) = (^)(°) + ( 93 2 )«, (6.1) 
where (c/? 2 )(°) is the vacuum expectation value for two parallel planes, and 



2\(1) 



i _ H oo 

X 2 



2 D -^b 



^^n? 1 siii 2 (n«)lif D /2-i l 



nu) 



n=l 



is induced by the presence of the plate at x = 0. In ((6.2(1 we use the notations 

u = 2irx/b, v = iry/b. 



Formula 1)6. 2J) can also be directly obtained by using the mode sum formula l)2T 
eigenfunctions for the boundary geometry under consideration: 



(fa 



2 sin(/cix) 
y /(2ir) D - 1 bu 



sin(-7rny/6) exp(ikrn — iuit), 



(6.2) 

(6.3) 
with the 

(6.4) 



where < k\ < oo, n = 1, 2, . . ., and u> = \Jk 2 + k\ + (Trn/b) 2 . 







y\ ! 








A 




N 


X 


\ 


S 



Figure 4: Geometry of two parallel plates with the interplate distance b perpendicularly inter- 
sected by the plate at x = 0. 

The representation similar to 1)6.1)1 takes place for the components of the energy-momentum 
tensor: 

(o|^|o) = (rn(°) + (^) {1) , (6.5) 

where the expectation values (Tj?)^ are induced by the intersecting plate at x = 0. They are 
related to the quantities investigated in previous section by formulae 

(T?)U=]hn(Ti) a , (^ 1 )W=-Umi(T 2 1 ) a , (6.6) 
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with lim corresponding to the limit a — > oo, (f)Q — > for fixed a — r and a<po- Using the formulae 
for (T^) a for the induced quantities one finds (no summation over i, i = 0, 3, . . . , D) 



7TX 



2 \ ~* -+1 



n=l 



(4e _l )(jD _l) 



nu 



sin (nv)K D / 2 (nu) 



+ 



2£ 



1 sin 2 (raw) 
- H — - 

2 D - 1 



(T[ 



{T a }i i) 



(T , 2) (i) 



vr(4£ - l)^ 1 -! 
tt(4£ - l)^ 1 - 



o OO 
2 



2^&f+ 2 n=1 



n^~ +1 cos(2nu )Er£)/2-i( n-u ); 
n=l 

oo 

n"2" +1 sin(2nu)i ; C£) / /2(^'u), 



(6.7) 
(6.8) 

(6.9) 



1-- 
7ra: 2 



n 2 



1 J (4£- 1)(£>- l) 



n=l 



sin 2 (nv)K D i 2 (nu) 



+ 



- + (4£- l)sin 2 (nw) 



K D /2-i(nu) 



(6.10) 



with the notations from (|6.3j) . Due to the presence of the MacDonald function in this for- 
mulae, the series are exponentially convergent. Note that in the alternative formulae for the 
energy-momentum tensor for minimally and conformally coupled scalar fields, given in the 
convergence of the series is power-law. 

For x > quantities H6.7j) - ()6.1l)|) are finite on the plates y = 0, b: 



y=0,b 



TTX 2 Ai 



n'^ x K D i 2 _ x {nu) 



(6.11) 



and (T^ ) 1 - ' \ y =o,b = 0. In particular, additional vacuum forces acting on the plates y = 0, b 
due to the presence of the plate at x = are determined by the component (T!^)^ 1 ' and are 
attractive. For large distances from the plate x = 0, by using the asymptotic formula for the 
modified Bessel function, we can see that the main contribution comes from n = 1 term and to 
the leading order we find 

(7f)« -iT^V", u»l, (6.12) 

with u defined by relation H6.3|) , In this limit the vacuum expectation values induced by the 
plate x = are exponentially suppressed. 



7 Conclusion 

We have investigated the Wightman function, vacuum expectation values of the field square 
and the energy-momentum tensor for a massless scalar field with a general curvature coupling 
parameter inside a wedge with a coaxial cylindrical boundary. We have assumed Dirichlet 
boundary conditions on the bounding surfaces. The generalization of the corresponding results 
for other boundary conditions is straightforward. The application of the generalized Abel-Plana 
summation formula for the series over the zeros of the Bessel function allowed to extract from 
the expectation values the parts due to the wedge without a cylindrical boundary and to present 
the additional parts induced by this boundary in terms of exponentially convergent integrals. 
The vacuum densities for the geometry of a wedge without a cylindrical boundary are considered 
in section |3J We have derived formulae for the renormalized vacuum expectation values of the 
field square and the energy- momentum tensor, formulae (|3.11|) . (|3.17j) . (|3.18|) . These formulae 
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can be further simplified for the important special case D = 3 [formulae (|3.15|) . (|3.25|) - (|3.28|) ]. 
For a conformally coupled scalar the energy-momentum tensor is diagonal and does not depend 
on the angular variable 4>. The corresponding vacuum forces acting on the wedge sides are 
attractive for 4>q < tt and are repulsive for (J)q > n. The previous investigations of the problem 
are restricted to this case. 

In section 0] we have investigated the additional expectation values for the field square and 
the energy-momentum tensor induced by the presence of the cylindrical surface. The field 
square is given by formula (|4.3|) and vanishes on the wedge sides <f> = <fi m for all points away 
from the cylindrical surface. The energy-momentum tensor induced by the cylindrical surface 
is non-diagonal and the corresponding components are determined by formula (|4.8j) . I|4.9|l . The 
off-diagonal component vanishes on the wedge sides. The additional vacuum forces acting on the 
wedge sides due to the presence of the cylindrical surface are determined by the ^-component of 
the corresponding stress and are attractive for all values <p$. On the wedge sides the correspond- 
ing vacuum stresses in the directions parallel to the wedge sides are isotropic and the energy 
density is negative for both minimally and conformally coupled scalars. The vacuum expectation 
values diverge on the cylindrical surface. For the points with \cf> — <p m \ 3> 1 — r/a the leading 
divergences are the same as those for a cylindrical surface without the wedge. As an illustration 
of the general results, we have plotted in figures [2] and El the components of the vacuum energy- 
momentum tensor induced by the cylindrical surface in the case of a conformally coupled D = 3 
scalar field for a wedge with the opening angle 0o = 7r/2. The Wightman function, vacuum 
expectation values of the field square and the energy-momentum tensor in the region outside 
the cylindrical shell are investigated in section [3 The formulae for these quantities differ from 
the corresponding formulae for the interior region by the interchange I qn (z) ^ K qn (z). For 
large distances from the cylindrical surface, r ^> a, the vacuum expectation values behave as 
(a/r) D ~ 1+2q for the field square and as (a/r) D+1+2q for the diagonal components of the energy- 
momentum tensor. As in the case of interior region, the vacuum forces acting on the wedge 
sides due to the presence of the cylindrical shell are attractive. In sectionH3we have considered a 
limiting case cfto — > 0, r, a — > oo, assuming that a — r and a<po are fixed. This limit corresponds to 
the geometry of two parallel plates perpendicularly intersected by the third plate (see figure 0J) 
and has been investigated previously in ^U] for special cases of minimally and conformally cou- 
pled scalar fields. For general values of the curvature coupling parameter the additional vacuum 
expectation values induced by the intersecting plate are given by formulae 1)6.2(1 . Q6.7JI - Q6.1UJ) . In 
particular, the corresponding energy-momentum tensor is diagonal on the parallel plates with 
the components Ijfi.llJ) and the additional vacuum forces due to the presence of the intersecting 
plate are attractive. Vacuum densities induced by this plate are exponentially suppressed for 
large distances. 

The generalization of the results obtained here for the Neumann, or more general Robin 
boundary conditions is straightforward. For instance, for the Neumann case in the expressions 
(|2.5|) of the eigenfunctions the function cos(qn(J)) stands instead of sm(qn(j>) and the quantum 
number n takes the values 0, 1, 2, . . .. In the case with a cylindrical boundary now the eigenvalues 
for 7a are zeros for the derivative of the Bessel function. The formula to sum the series over these 
zeros can be taken from |lfi| . Now, in formulae (|4.3[) . (|4.8jl . 1)4. 9|) for the vacuum expectation 
values, instead of the ratio K qn {z) / I qn {z) the ratio of derivatives, K' qn (z) / I' qn {z) , will stand. 
Note that in this paper we have considered boundary induced vacuum densities which are finite 
away from the boundaries. As it has been mentioned in Ref. [2B1; the same results will be 
obtained in the model where instead of externally imposed boundary condition the fluctuating 
field is coupled to a smooth background potential that implements the boundary condition in a 
certain limit |29j . 
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